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, Abstract. We give a progress report of our lattice calculation of direct and indirect CP violation in 

kaon decays, parametrized as e'/ e and Bk, which require non-perturbative calculation of the matrix 
elements of the Standard Model effective Hamiltonian. 

m : 

\D : INTRODUCTION 

■ In the investigations of the Standard Model, a very important issue is the theoretical 

treatment of K — > %% decay to the accuracy such that comparison with the experimental 
results is possible. In particular, the ratio of direct and indirect CP violation, e'/e, has 
been determined experimentally |]]] in recent years and theoretical calculation is desired 
to test Kobayashi-Maskawa theory. In the theoretical calculation, numerical simulation 
of lattice QCD is the most systematic method to estimate the non-perturbative effect of 
QCD which is the main source of the error. Using the operator product expansion, the 

interaction in this decay is written as H w = GFV /J^ ud W; (jU ) Qi , where the coefficients 

Wj contain the effects of the energy scales higher than the matching point \x and can be 
obtained perturbatively |y]. Non-perturbative QCD effects will appear in K — > kk matrix 
elements of the local operators {nit \Qi\K), which should be calculated on the lattice. A 
couple of years ago, CP-PACS and RBC Collaboration flUl] calculated all of the matrix 
v-j ■ elements using the domain-wall fermion formalism o,|g] to realize the chiral symmetry 

required in this calculation and reported small and negative values of e'/e in conflict 
with the experimental result. Another work using staggered fermion has obtained a larger 
negative value [7]. In these works, however, there are several uncontrolled systematic 
errors such as 1) the effect of the small, but non-zero, chiral symmetry breaking, 2) the 
effect of finite lattice spacing, 3) the effect of the perturbative treatment of the charmed 
quark in the matrix elements, 4) quenching effect, and 5) K — > kk matrix elements are 
obtained from K — > % and K — > (vacuum) by using lowest order chiral perurbation 
theory & 

In order to examine all of these systematic errors except the fifth one, we are perform- 
ing two types of numerical simulation with domain-wall fermion and the DBW2 gluonic 
action [9] to improve the chiral symmetry on the lattice. "Numerical Simulation I" is the 
quenched calculation including directly the effect of the charm quark on the lattice. The 
degree of chiral symmetry breaking is decreased by a factor 1/10 compared with the 
previous work of RBC Collaboration. In addtion, we are generating gauge configura- 
tions with Nf = 2 dynamical quarks HlOfl in "Numerical Simulation II." In the rest of this 



article, we present the contents of these numerical simulations and report preliminary 
results of the matrix elements which numerically dominate e' / e and kaon B-parameter 
Bk- 



NUMERICAL SIMULATION I 

We are generating gauge configurations on a relatively fine 24 3 x 48 lattice with the 
scale cT x = 2.86(9). The residual quark mass m res which measures the chiral symmetry 
breaking is as small as < 0.3 MeV. Since quark mass mja is introduced as a parameter 
of the boundary condition in the fifth dimension in domain-wall QCD, the localization 
of chiral modes on domain-walls in the fifth dimension tends to fail for a heavy quark 
mass nif. However, our small lattice spacing made the value of m c a acceptable as a 
domain- wall fermion: m c a ~ 0.45. We found that, around this value, the behavior of 
wave function in the fifth dimension is qualitatively same as the case of much smaller 
quark mass. 

At the lowest order of chiral perturbation theory, K — > %% matrix elements are in 
proportion to K — > % matrix elements calculated on the lattice. For i = 1 - 6, 9, 10, these 
matrix elements are related as, 
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in particular, for AI = 1/2, or / = 0, subtraction of a lower dimension operator is needed: 
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For i = 7, 8, we have the simpler relation 
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In particular, and have the largest contribution to e'/e, numerically. In FIG.d 
results of K — > it matrix elements of these operators are plotted. In particular, the left 
panel, which is the example with m c a = 0.40, shows that there is a severe cancellation in 

the subtraction in © for . Since the slope of the subtracted matrix elements have the 
error of ~ 200%, we cannot quote a result of K — >• %% matrix element with the current 
statistics. And its depencence on m c a is not visible, so far. 
Lattice value of kaon B parameter which is defined by 
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FIGURE 1. K -> % matrix element of (left) and ) (right) as a function of m fa from Numerical 
Simulation I. In the left panel, data for matrix elements before (circle) and after (diamond) the subtraction 
and the subtraction term —a,(,(n |2sub|^) (square) are plotted from 50 statistics. Linear extrapolation was 
used for all of plots. 
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FIGURE 2. Lattice value of Bk as a function of mta from 77 configurations. 



is plotted in FIG. El as a function of m^a. In this figure, the fit function used is Bk = 
£o[l +Cmfaln(mfa)] + ^imfa with C taken from analytic result 111 ill . The physical result 
for Bk can be obtained at nif = m s /2 (the filled symbol). To obtain the physical value, 
we are now calculating the Z factor for Bk by non-perturbative renormalization proposed 
in lfl2ll and the preliminary result is roughly consistent with the previous works lfl3L Ell . 



NUMERICAL SIMULATION II 



Since the dynamical simulation demands much more resources than a quenched one, 
dynamical domain- wall QCD has been explored by only our collaboration ill Oil , so far. 



5 




00 




0.4 
0.3 


• unsubtracted si-^ 

• subtraction term 

• subtracted 


-0.01 
-0.02" 




0.2 




-0.03 




0.1 

o.o = 


~ss — • • ^ 


-0.04 
-0.05 
-0.06 




0.1 







I , 1 , 1 , 1 , 1 , 1 , 1 _ 06 I , 1 , 1 , 1 , 1 , 1 , 1 

0.01 0.02 0.03 0.04 0.05 0.06 0.01 0.02 0.03 0.04 0.05 0.06 



m f a m,a 

FIGURE 3. Same as FIG. ^but from Numerical Simulation II with m sea a = 0.03. 72 configurations 
were used. 



In this calculation, we generated three kinds of gauge configuration on a 16 x 32 lattice 
with the mass of the sea quark (u or d quark) being m sea a = 0.02,0.03 and 0.04. For 
each series of configurations, K — > % and K — > matrix elements are calculated in the 
same way as Simulation I with the five valence quark masses m va i = 0.01-0.05, and 
basic parameters cT l 1.8 GeV and m res « 3 MeV are obtained. FIG.|3]shows the same 
matrix elements as in FIG. [T]as an example of the case of m sea = 0.03. Although the 
signal seems to be reasonable, we need much more statistics to take correct chiral limit 
m sea a = m w& \ — > using three data points with m sea a = m va i = 0.02, 0.03 and 0.04. B K at 
the physical point such that m sea _= m u u and m s ~ 120 MeV will be also obtained after 
a careful treatment of our data II14II and non-perturbative renormalization which is now 
under calculation. 

We thank RIKEN, BNL and the U.S. DOE for providing the facilities essential for the 
completion of this work. 
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